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On the velocity dispersion of young star clusters: super- virial or 
binaries? 
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ABSTRACT 

Many young extra-galactic clusters have a measured velocity dispersion that is too high for 
the mass derived from their age and total luminosity, which has led to the suggestion that they 
are not in virial equilibrium. Most of these clusters are confined to a narrow age range centred 
around lOMyr because of observational constraints. At this age the cluster light is dominated 
by luminous evolved stars, such as red supergiants, with initial masses of ~13— 22 M Q for 
which (primordial) binarity is high. In this study we investigate to what extent the observed 
excess velocity dispersion is the result of the orbital motions of binaries. We demonstrate that 
estimates for the dynamical mass of young star clusters, derived from the observed velocity 
dispersion, exceed the photometric mass by up-to a factor of 10 and are consistent with a 
constant offset in the square of the velocity dispersion. This can be reproduced by models of 
virialised star clusters hosting a massive star population of which ^25% is in binaries, with 
typical mass ratios of ^0.6 and periods of ^1000 days. We conclude that binaries play a 
pivotal role in deriving the dynamical masses of young (^lOMyr) moderately massive and 
compact (< 10 5 M Q ; > 1 pc) star clusters. 

Key words: globular clusters: general - open clusters and associations: general - galaxies: 
star clusters - binaries: general - binaries: spectroscopic - supergiants 



1 INTRODUCTION 

Young massive clusters have received considerable attention 
in the last decade because they trace star formation (e.; 



IWhitmore & Schweizedl 19951 ; iMiller et alj|l99l IZepf et al.ll 1999 
Advances in observations enabled us to resolve such star clusters 
up to ~ 20 Mpc, allowing determinati on of th eir fundamental pa- 
rameters, such as mass and radius (e. g. lLarser]|2004) . 

The mass of a resolved star cluster can be determined in two 
ways: one of them by converting the observed luminosity, age and 
distance directly to mass via the age dependent mass-to-light ra- 
tio (M / L) taken from a single stellar population (SSP) model. We 
refer to the resulting mass as the photometric mass Afphot. This 
method requires an estimate of the cluster age, which again requires 
estimates for the metallicity and the stellar initial mass function 
(IMF). 

An independent mass estimate is based on the virial theo- 
rem and this mass is generally referred to as the dynamical mass 
dSpitzedl 19871) : 



dispersion in the cluster, r e tt is the effective (half-light) radiu^] and 
rj ~ 9.75 is a constant that depends slightly on the density profile. 

Equation l[T) is valid for a cluster in virial equilibrium consist- 
ing of single stars. Since in this study we consider possible differ- 
ence between crd yn and the observed velocity dispersion, rj bs, we 
will refer to the empirically derived dynamical mass, i.e. based on 



fobs 



, as AC 



dyn " 

A comparison between Af p hot and ACj™ serves as a check for 
the range of assumptions on which both mass estimates are based. 
An inconsistency between M p hoi and Mj 5 ™ can be attributed to vari- 
ations in the IMF, on which Af p hot is in part based, or to a lack 
of virial equilibrium, on which is based. For many young 

(~10Myr) star clusters Afff > M phot , with M H ° v b n s up to ~10 
times larger than Afphot ( e -g- Bastian et al.ll2006l hereafter B06), 
suggesting that these objects are super-virial. For older clusters 
Gt, 100 Myr) there is good agreement between ACS and Af p hot (e.g. 
iLarsen et alj|2004 B06). 

The alleged super-virial state of some young clusters has 
been attributed to the impulsive expulsion of residual gas from 



(1) 



Here G is the gravitational constant, <Td yn is the line of sight velocity 



1 Here we assume that the half-light radius is the same as the h alf-mass ra- 
dius, which is not t he case whe n the cluster is mass segregated (Fl eck et all 
l2006tlGaburov & Gielesll2008l) . 
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the parent molecular clou d in which the star cluster formed (e.g. 
iGoodwin & Bastianl2006h . Such early outgassing, driven by stellar 
winds of massive stars or supernovae, causes the stellar velocities to 
be high compared to the binding energy of the stars. This argument 
has been used to m otivate infant mortality of young star clusters 
dLada & Ladall2003l) . 

However, the gas expulsion theory has difficulties in explain- 
ing the super- virial velocity for the lOMyr old clusters presented 
in B06. The arguments are as follows: The time needed to com- 
pletely dissolve, or to find a new virial equilibrium after impulsive 
gas expulsion is about 20 crossing times, t a , where t a oc p^ 1 ^ 2 
and ph is the density within the half-mas s radius (see for exam- 
ple Fig. 8 in Baumg ardt & Kroupal 12007). Hence, to be able to 
'catch' an unbound or expanding cluster at 10 Myr, t a should be 
> 1 Myr. This corresponds to a half-mass density of stars and 
gas of ph < 300 Mq pc -3 . Clusters with shorter t a (higher den- 
sity) have expanded into the field, or found a new equilibrium, 
a few Myrs after gas expulsion and are not observable as super- 
virial clusters at lOMyo The density in the embedded phase of 
the clusters under discussion is unknown, but can be roughly esti- 
mated using their current densities. The present day densities are 
p h ~ 10 3±1 M Q pc -3 (Table [3}. The densities in the embedded 
phase were at least a factor 1/e 4 higher, where e is the star for- 
mation efficiency. This because the mass of the embedded cluster 
has reduced by a factor e and the cluster has expanded at least by 
a factor of 1/e as a response to it, contributing a factor 1/e 3 to the 
reduction of pi,. The 1/e expansion holds for adiabatic mass loss, 
for impulsiv e mass loss and e < 0.9 the cluster expands much more 
jHillsl 1980h . So for the clusters at 10 Myr the estimated densities in 
the embedded phase are much too high to still have features of gas 
expulsion detectable in their velocity dispersion at 10 Myr. These 
arguments suggest that deviations from virial equilibrium are not 
a plausible explanation and an alternative explanation for the high 
o" bs values is needed. 

The existence of binary stars is generally ignored in the es- 
timates for M^l, even though their internal velocities can lead to 
an over estimation of Md yn (Kouwenhoven & de Grijs 2008, here- 
after K08). K08 studied this phenomenon in virialised star clusters 
with a 100% binary fraction and a range of (Jd yn . They subsequently 
derive M^l by 'measuring' <7 b s and applying equation 10. They 
found that the presence of binaries can lead to an overestimation 
of Aidyn by a factor of ~2 for clusters with <7dyn — lkms -1 . For 
clusters with <Jd yn — lOkms -1 they found only a 5% increase in 
M^yn due to binaries. They therefore c oncluded that binari es are 
not important for massive/dense clusters. fMengel et alJd2008l here- 
after M08) found A/ d ° y b n s /M p hot — 10 for some of the star clusters in 
the Antennae galaxies (NGC 4038/4039) and NGC 1487 and since 
these clusters have velocity dispersions of 10 — 20 km s _1 they sub- 
sequently concluded that binaries are not important and that these 
star clusters are super-virial and dissolving quickly. 

Here we revisit the effect of binaries on iWJ™ /M p hoi and we 
focus on ~ 10 Myr old star clusters. This is motivated by our desire 
to incorporate the effect of the steep increase of the stellar lumi- 
nosity with increasing stellar mass, which is in particular important 
for young clusters, an aspect not considered by K08. Of the two 
approaches presented by K08, one focused on solar-type stars, the 
other uses a Kroupa IMF for the primary stars. They give equal 



2 The models of IGoodwin & BastiarJ 120061) start with a density of ~ 
60 Mq pc — 3 (t C r & 2.5 Myr) in the embedded phase, and this is why they 
find that the effects of gas expulsion are observable for 25 Myr. 



weight to each binary in their computed velocity dispersion. In ad- 
dition K08 do not consider stars more massive than 20 Mq . This 
approach may be appropriate for studies of intermediate age and 
old open star clusters, but it is less suitable for young star clusters. 

At an age of ~ 10 Myr the cluster light is dominated by 
the most massive (> 15 Mq) stars for which binarity is high 
and ignoring them can lead to misinterpretat ions of observa- 
tions of various astrophysical processes (e.g. IVanbeveren et al.l 
1998). Massive binaries have a larger effect on a bs than low- 
mass binaries due to their higher orbital velocities, but also due 
to the more common shor t -periods and comparable masses (e.g. 
Duquennov & Mayor 1991; lPqrtegies Zwart et al . 2002; Sa na et al.l 
2008, 2009; Mas on et alj|2009ri . Incorporating the massive stars in 
our calculation has two important effects, both of which amplify the 
effect of binarity on <r b s with respect to the results of K08: massive 
stars dominate the cluster light and their higher masses and (intrin- 
sically) different binary properties give rise to a larger a b s . 

In this paper we quantify the effect of the presence of (mas- 
sive) binaries on M^/Mphot and we use this ratio is a proxy of the 
excess dispersion. In Section[2]we discuss the properties of the bi- 
nary population that is expected in young (~ 10 Myr) clusters and 
we present a simple model for the additional velocity dispersion 
due to such binaries. In Section [3] we summarize existing obser- 
vational results to confront our model with. Our conclusions are 
discussed in Section|4] All the specific acronyms used in this study 
and their definitions are given in TableQ] 



2 THE VELOCITY DISPERSION DUE TO BINARITY 
2.1 The importance of massive binaries 

The young clusters with measured and M p hot have a rather 
narrow range in ages of ~8 — 13 Myr. This is mainly because of the 
onset of red supergiants in this age range making clusters brighter 
and easier to detect and study in detail. Stars in a stellar popu- 
lation with an age of 10 Myr have initial mass of 13 — 22 Mq, 
corresponding to masses of 13 — 16 Mq at an age of 10 Myr 
dLeieune & Schaerer 2001). If we would consider a small spread 
around 10 Myr the quoted mass range would be slightly larger, but 
for simplicity will continue with the assumption of a constant age 
of 10 Myr. Those massive stars appear to have high primordial mul- 
tiplicity with a spectroscopic binary fraction of ~50% or more (i.e. 
/ > 0.5, lMason et al.ll2009l ; lBo"sch et al.ll2009h . 

Most of the measurements we discuss in Section [3] are done 
in the near infrared. At these wavelengths red supergiants (RSGs) 
dominate the observed light and therewith the measured (Jobs- For 
the studies done in the optical wavelength the blue supergaints 
(BSGs) are more dominant. We here refer to the population of lu- 
minous evolved stars as supergiants (SGs) and use the subscripts 
SG to denote parameters that apply to these stars. 

Since the SGs outshine the main sequence stars it is important 
to establish the binary fraction among them. This estimate is com- 
plicated by the internal evolution of binary stars affecting especially 
the RSG phase and hence the actual population of SGs present at an 
age of 10 Myr. In particular, short-period binaries are likely to ex- 
perience a common envelope evolution (CEE) and/or Roche-lobe 
overflow (RLOF) which causes the binary components to follow a 
different evolution compared to single stars of similar initial mass, 
and may prevent the RSG stage altogether. 

lEldridge etakl d2008l hereafter E08) find that these effects re- 
duce the average duration of the RSG phase by a factor of two or 
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three. They find this for a population of binaries with a flat distri- 
bution in log P/d between —0.15 and 4.5 and a flat distribution 
of q between 0.1 and 0.9, where P and q are the orbital period 
and the ratio of the secondary mass over the primary mass, respec- 
tively. The short period binaries with high mass ratios are most 
affected by interactions through RLOF and CEE. For our simple 
model we assume a minimum period, P cl ;t, and as an approxima- 
tion of the shortened evolved phase of primaries in tight binaries 
we remove the binaries with P < P cr j t . The fraction of binaries we 
remove should roughly match the fractional reduction of the aver- 
age life-time of the RSGs (factor of 2 — 3). This constraint is met 
for Pcrit = 500 d since 63% of the binaries in the E08 population 
haveP < 500 d for mi = 15 M . 

RSGs at lOMyr have a maximum radius of ~900Rq. For 
mi = 15 Mq and q = 0.6 this corresponds roughly to the separa- 
tion of a binary with P cr i t . For P = 20 00 d the Roche- lobe radius is 
around 900 R© (using the formula of Eggl etonll 19830 and binaries 
with longer orbital periods will follow an evolutionary path sim- 
ilar to single stars (08). So in our model we remove all binaries 
with P < 500 d and assume that binaries with P > 500 d ex- 
perience a SG phase unaffected by binary evolution, even though 
it is expected that the RSG phase of primaries in binaries with 
500 < P/d < 2000 is affected by the companion. In reality it 
will not be such a step function, since most SGs do contribute at 
some stage in their evolution to the integrated light. But under our 
assumptions we reduce the number of binaries roughly by the same 
fraction as what was found for the fractional reduction of the av- 
erage RSG phase in the model of E08. By removing all binaries 
with P < 500 d we are probably making a conservative approach 
since we bias our binary population to longer periods. In reality 
these binaries can continue to contribute to the velocity dispersion. 
This because the primary does not necessarily becomes dark af- 
ter its shortened RSG phase, and if it does, the secondary can still 
contribute to the velocity dispersion (E08). 

The relevant parameter for studying the binaries that con- 
tribute to the velocity dispersion is the fraction of binaries among 
SGs, which we identify with /so - Using N p for the number of stars 
with initial masses in the range 13 — 22 Mo and the fraction of 
binaries with an orbital period P > P c ,i t as g, then the number of 
stars in binaries unaffected by interaction is gfN p and the number 
of SGs that is removed is (1 — g)fN p . The total number of remain- 



ing SGs, i.e. single and in binaries, is (1 
can write for fsa 



fsG 



gf 



(!-/) + • 



f)N t + gfN t . So we 



(2) 



In equation lO we have neglected the possibility that secondary 
stars contribute to the SG population, thus slightly underestimating 
fsa- If all stars are in binaries (/ = 1) then fsa = f for all values of 
g. For the remainder of our analysis we adopt a more conservative 
value of / = 0.6 in our parametric model (Section[Z2j and a range 
0.3 < / < 0.9 for the Monte Carlo simulations in Section [2~3l 

The orbital periods of early-type spectroscopic binaries range 
from a couple of days to about 10 years. Adopting an Opik's law 
in the interval 0.3 < log P/d < 3.5 and a period threshold 
logP cl it/d = 2.7 we find g — 0.25 (i.e. we remove 75% of the 
binaries), which via equation $2% results in fsa — 0.25. For the 
representative period we use P = 10 s d, which is approximately 
the mean of the periods above P cl -i t when assuming a flat distribu- 
tion in log P. 

The distribution of mass ratios for high-mass stars appears to 
be flat between q ~ 0.2 (the typical detection limit for SB2 sys- 



Table 1. Overview of the specific acronyms used in this study. 

Acronym Description 

/ primordial binary fraction of massive stars (13 — 22 Mq) 

g fraction of primordial binaries unaffected by interaction 

/sg effective binary fraction among SGs at lOMyr 

mi mass of the primary star 

N p number of stars with initial masses in the range 13 — 22 Mq 

q ratio of the secondary mass over the primary mass 

P orbital period 

P cr i, minimum period for binaries to be unaffected by interaction 

i? or b orbital velocity of the primary star 

i?id line of sight velocity of the primary star 

r e ff cluster half-light radius in projection 

crdyn ID dynamical velocity dispersion of cluster members 

tr bs empirically determined ID velocity dispersion 

crbin ID velocity dispersion due to binary orbital motions 

dyn dynamical cluster mass based on a^ yn 

empirically determined dynamical mass based on cr3 

A^phot photometric cluster mass 



terns) and q = 1 (e.g. lSana & Le Bouquml2 009). We adopt q = 0.6 
as a typical value for the mass ratio. 

Our adopted values of the parameters that control the SG bi- 
nary population at lOMyr are summarised in Table[2] These values 
serve as input for the model presented in the next section. 



2.2 A parametric model for the velocity contribution of 
binaries 

To quantify the importance of binaries on cr bs we model their ob- 
servational characteristics. Since the dynamical velocities of the 
cluster members (stars and centres of mass of binaries) and the 
orbital velocities of the binary members are uncorrelated, we can 
write <7o bs = a| yn + <r bin . Here we derive a simple expression for 
the contribution to a^ hs of the orbital motions of binaries, <7 bin . 

Since the secondary is generally much fainter than the primary 
we ignore its contribution to the light and focus only on the primary 
star. Its orbital velocity, w rb, can be expressed in terms of q, mi and 
P using Kepler's third law: 



Vovb = q 



1 + 9 



2/3 / si \ 1/3 

2P 



(3) 



The contribution to the line of sight velocity, vid, depends on 
the inclination, i, of the orbital plane and the phase, 9, in which 
the binary is observed. We first assume a population of binaries 
with the same q, mi and P and random orientations of the orbital 
planes and (un-correlated) random orbital phases. This results in 
flat distributions of —1 $J cos(i) ^ +1 and ^ 9 ^ 2ir. For 
each individual binary vm = v m b sin(i) cos(f?) so the distribution 
of vid values is the joint probability density function of sin(i) and 
cos(f?) multiplied by n olb , which is flat between — v rb and +foib- 



The variance of this distribution is <r bir 



b /3. In reality there 



will be a spread in the binary parameters which will make the line 
of sight velocity distribution peaked, with a similar variance. We 
continue with the assumption of a population of identical binaries 
to be able to analytically express our result in the binary param- 
eters. In Section [231 we validate this assumption and quantify the 
expected spread using Monte Carlo simulations. 

Taking into account that only a fraction fsa (Section [2.1 1 ) of 
the stars that contribute to the cluster light is part of a binary re- 
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Table 2. Adopted values for the parameters of SG binaries at lOMyr. 





Reference 




Range 


parameter 




min 


max 


distribution 


/ 


0.6 


0.3 


0.9 


flat 


9 


0.25 








fsG 


0.25 








mi/M 


15 


13 


16 


Salpeter 


1 


0.6 


0.2 


1.0 


flat 


logP/d 


3.0 


0.3 


3.5 


flat 


logP crit /d 


2.7 


2.7 


2.7 





duces (j 2 in by a factor /sg- The dependence of <r 2 in on the binary 



parameters can then be expressed as 

,3/2x4/3 



2 

"bin — 



JSG 



2q* 



TvGmi 
2P 



2/3 



(4) 



For the reference values (Table [2]l we find that ab m — 
6.6kms _1 , which is equal to crdyn for a (virialised) cluster with 
M = 10 s Mq and r c ff = lpc (equation QJ. So for such clusters 
and these binary parameters overestimates the true mass M 
by a factor of two because of binaries. We use these scaling values 
to write a more general expressing for the ratio 



_2 

°bin 



'dyu 



fsG 

0.25 



q \i 

0.6/ 



/ mi 
\15Mq 



' / M/rgt 
VlO 5 M pc- 



where we have approximated the term [2g 3//2 /(l + g)] 4 ^ 3 
equation © by g 3 ^ 2 . Equation {5} is accurate to within 8' 

q Z 0-2. 



, (5) 

from 
h for 



In the next section we will use the ratio M^l/M as a measure 



of the excess dispersion, which we can write as 



M 



2 2 
°"dyn + ""bin 



'dyn 



1 + 



M/reff 



(6) 



(7) 



lO 5 M pc" 

where in the last step we have used the reference values of Table [2] 
such that the binary part of equation lf5J equals 1 . 

For M/reff < 10 5 Mq pc _1 binaries dominate the measured 
velocities and therefore M^l/M oc (M/r e n)~ l (for a constant 
cry,,)- For higher values of M/r e « the presence of binaries has little 



effect on the estimated mass and M^l/M : 



1. 



2.3 A Monte Carlo validation 

Up to this point we have assumed populations of equal binaries 
giving a flat distribution of ?;id values and a fixed value for own for 
each cluster. First we verify the assumption that the shape of the 
v\o distribution resembles a Gaussian when a range of binary pa- 
rameters is assumed (Section [2.3.1l l. Then we quantify the expected 
spread in <r 2 in values when comparing clusters (Section[2372]l. 



2.3.1 The velocity dispersion of a binary population 

We generate two populations of 10 4 binaries, i.e. no single stars, 
to study the shape of the velocity dispersion of their orbital mo- 
tions. For one population we give all primaries a w or b based on equa- 
tion <[3j and the reference values from Table [2] The values of vw 



o.o4 r 



0.03 



0.02 - 



0.01 - 



Reference values 
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0.00 




-30 -20 -10 

1> 1D [km s~ 
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20 



30 



Figure 1. Monte Carlo simulations of the probability density functions of 
■Uirj values of binary populations based on the reference values (i.e. all bina- 
ries identical, dotted histogram) and using a spread in the binary parameters 
(dashed histogram). Both simulations are based on the values for mi, q and 
logP quoted in Table [2] and both consist of 10 4 binaries. The Gaussian 
curve shown with a full line is the approximation based on the reference 
values from Tableff] but using / = 1 (cr bin = u or b/\/3 13kms -1 ). 
This approximation nicely describes the more realistic simulation based on 
a range of values (dashed histogram). 



are acquired by multiplying v 0I \, for each binary by a random num- 
ber between —1 and +1 (Section [2.2b . The resulting distribution 
and the Gaussian approximation (obm = v OT b/V3, SectionffTJJ are 
shown as a dotted histogram and a full line, respectively, in Fig.fTJ 
For the second population we randomly draw values for the masses, 
mass ratios and periods from the distributions described in TablefJ] 
The binaries with P < P cl j t are taken out of the sample. With equa- 
tion {3]l we then calculate v or b for each remaining binary and bid 
is again acquired by multiplying v m b by a random umber between 
— 1 and +1. The resulting distribution is shown as a dashed his- 
togram in Fig. [TJ Two things can be seen from this figure: 1.) the 
width of the more realistic distribution (i.e. using a range in binary 
parameters) is well approximated by our simple model and 2.) this 
distribution is close to Gaussian. This last point is important since 
we have assumed in Section [2~2l that we can quadratically add Obm 
to Ody,, to get the total velocity dispersion. 



2.3.2 The expected dispersion in the binary dispersion 

Here we quantify the spread in <T 2 in > i.e. the dispersion in the addi- 
tional velocity dispersion squared, when comparing different real- 
isations of binary populations, caused by the fact that the number 
of binaries is small and that there is a spread in the binary fraction 
(TableO. 

We generate 1000 massive star populations, each consisting 
of 200 SGs (an ap proximate number for a cluster of mass 10 5 Mq, 
lLarsen et al.l 120081) . For each population, we randomly sample a 
value for / and thus have 200 x / binaries. The ?;id values of the 
binaries are calculated in the same way as in Section [2. 3 . 1 1 using 
the ranges from Table[2] For each population the variance of the 1- 
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Figure 2. The probability density function of a2_ values following from 
the Monte Carlo experiment described in Section [23] for the adopted dis- 
tributions in /, mi, q and log P (Table|2j. The solid line near the peak of 
the distribution indicates the value of a?, derived with equation (4} using 
the reference values from Tableff] The dashed lines indicate factors of two 
of variation. This corresponds approximately to the one and two standard 
deviations (horizontal arrows) of the log-normal approximation. 



dimensional velocity distribution (o" bjn ) of the remaining SGs in the 
sample (single and binary) is then calculated. The resulting proba- 
bility density function (PDF) of the oy n values is shown in Fig. [2] 
The reference value of <r bin is indicated with a vertical solid line and 
is very close to the mode of the distribution. When approximating 
PDF(cr b 2 in ) by a log-normal we find a standard deviation of ~0.7 
corresponding to a factor of ~2 relative to the mode. 

We will now compare the model to empirical determinations 
of M$*/M. Since we do not know the real mass M we use M p hot 
as a proxy. 



3 COMPARISON WITH OBSERVATIONS 

We illustrate the effect of the presence of binaries by comparing 
the results of our model from the previous section to the empirical 
ratio Mj™/JVfphot f° r a number of clusters. The cluster masses M<j5 
and A/phot follow from literature values for magnitude, age, cr bs and 
r e ff. We subsequently re-derive M p hot and to obtain a homo- 
geneous sample, which is important because the literature values 
are derived by a number of groups using a variety of SSP models 
to derive Mphot and apply different (small) corre ctions to the value 
of rj (equation QJ because of mass segregation ( iFleck et alj|2006h 
and variations in the density profiles (for those clusters for which 
measurements of their surface brightness profile are available). All 
cluster parameters and references to the relevant literature are given 
inTable[3l 

We use the iBruzual & Charlotl d2003h SSP models with a 
Chabrier IMF and solar metallicity to derive M P h t- For each clus- 
ter the age dependent M/ L is found from the observed cluster age. 
Combining M/L with the absolute magnitudes (My for 7 clus- 
ters, and Mk for the rest) we determine Mphot. The quoted up- 
per and lower limits in M p hot are calculated through the uncertain- 



Figure 3. The ratio M^/M as a function of the square of the velocity 
dispersion expressed in terms of the observables Mphot A*eff (equation [TJ. 
The solid curve is calculated assuming the fixed binary parameters adopted 
in Section [2TT1 which are: /sg = 0.25, q = 0.6, m\ = 15 Mq and 
P = 1000 days. The dotted curves are calculated by varying a?, with 
factors of two for each subsequent curve. The symbols with error bars are 
the observed values for these parameters from TablefJ] The horizontal line 
is plotted to guide the eye. 



ties in log(age/yr). We use equation (TJ to determine Mj™, with 
rj = 9.75, and the uncertainty is calculated using the uncertainties 
in (Jobs and r^a by adopting standard error propagation. The calcu- 
lated values for Mphot and Mj™ are presented in Table [3] 

We now assume that our choice for the IMF and the metallic- 
ity is representative for all clusters and that variations in r\ due to 
mass segregation and the density profile are negligible. Under these 
assumptions Mphot reflects the true mass M and subsequently cr| yn 
scales with Mphot A"eff (equationQJ. However, M p h« is also affected 
by binarity since the shortened RSG phase in short pe riod binaries 
reduc es the integrated luminosity (recently noted by Davi es et al.l 
2009). This effect reduces the fraction of bright stars visible at 
lOMyr by a factor of (1 - g)f w 0.5 (SectionrXTl 

In Fig. [3] we present the data. The trend that clusters with a 
small Mphot/?"eff tend to have high Mj™ /Mphot, and which drops 
with increasing Mphot/Veff is well reproduced by a population of 
binaries among the most massive stars. The dispersion in the obser- 
vations around the mean value for our model (equation (7), solid 
curve in Fig. [3) roughly corresponds to the spread following from 
our Monte Carlo experiment (dashed lines) when allowing a spread 
in the binary parameters, rather than fixed values. 

The linear Pearson correlation coefficient 

( Rodg ers & Nicewanded Il988l) for the logarithmic values of 
the data presented in Fig. [3] is s = —0.71 with a significance 
level (s.l.) of ~2x 10 -4 , which indicates that the observed trend is 
statistical significant^. The downward trend in Fig. [3] suggest that 



3 The coefficient s can have a value between —1 and +1, where — 1(+1) 
indicates a linear relation between the observed variables with nega- 
tive(positive) slope. A value of s = indicates a lack of correlation. 
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Table 3. Overview of th e observational data use d. Values for the a bsol ute magnitudes, lo g(age/yr), <T obs and r e ff wer e take from (1) Ss tlin et alj J2007ri . 
f2) [Moll et all 120071): OllLarsen et al] bOQgl): f4l|Smifh et alj 120061) : r5) lMcCradv & Gra ham (2007); (6l lBastian eilul feOOd B0 6) and references th erein; 
m lMcCradv et alj J2003I) ; (8) lMengel et alj bOOa M08). The My value of NGC 6946-1447 was taken from the update given bv lLarsen etat] J20Q6I) . The 
values for MSjj and iWphot were re-derived in this study, see Section[3]for details. 
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a% bs equals cr| yn plus a constant. This is what follows if all clusters 
are virialised and host a similar binary population (equation[|>]l. 



4 CONCLUSIONS AND DISCUSSION 

Several studies have found from spectroscopic analyses that for 
many young (~10Myr) star clusters the measured velocity disper- 
sion is too high for the mass derived from their total luminosities 
and their ages. This has led several authors to conclude that these 
clusters are super-virial and thus dissolving. However, the conver- 
sion from velocity dispersion to mass (equation Q} does not con- 
sider the additional velocities of binaries. K08 considered this ef- 
fect, but concluded that binaries are only important for clusters with 
low intrinsic velocity dispersion (~lkms _1 ), i.e. lower than the 
aforementioned clusters. K08 ignored the mass dependent mass-to- 
light ratio of stars and the intrinsically different binary properties of 
massive stars. In this study we show that taking these aspects into 
account makes the contribution of binarity to the dynamical mass 
estimates, iWj™, of clusters in this age range non-negligable. 

We present a simple analytical model that gives the 1- 
dimensional velocity dispersion of a virialised star cluster hosting 
a binary population. The model is complementary to the classical 
virial relation for clusters consisting of single stars (equation [T}. 
The result is presented as a single equation that needs as input the 
(typical) binary fraction, mass ratio, primary mass and orbital pe- 
riod of the binary population and the mass and radius of the star 
cluster. This relation can be used to easily estimate the effect of 
binaries based on different parameters for the binary population 
and/or cluster. The model presented here serves as a starting point 
for more realistic approaches using binary population synthesis 



models (e.g Eldri dge & Stanwavll2009h . Tentative confirmation of 
our results comes from the velocity dispersion of the binary popu- 
lation discussed in E08: ~12 km s -1 at an age of lOMyr (Eldridge 
2009, priv. comm.), which is close to what we find for the reference 
values discussed in Section[2](see Fig.QJ. 

For 24 clusters we derive the ratio of M^* over the photomet- 
ric mass, Mphot! an d show that it decreases with increasing cluster 
velocity dispersion. This is also what follows from the model and 
most of the empirically determined M^™ /Mphot ratios can be ex- 
plained by binaries using a conservative binary fraction of 25%, a 
mass ratio of 0.6 and an orbital period of a 1000 days. When allow- 
ing a spread in the binary parameters, almost all clusters are within 
2-standard deviation of the model results. 

The fact that and A/ p hot generally agree for older (> 

100 Myr) clusters is consistent with this binary scenario. In older 
clusters, we indeed expect a lower velocity contribution of binaries. 
The primary star will be of a later spectral type, thus mi is lower. At 
100 Myr the most luminous stars are roughly 5Mq. Equation I0 
shows that when mi is a factor of 3 lower, a bm is a factor of ~2 
lower, keeping all other parameters fixed. Also, typical periods are 
longer. Duquennoy & MayoJ ( 1 199 ll) find that the median period of 
solar type stars is 180yr. From equation (O we can see that the ef- 
fect of such binaries on a{; m is about a factor of ~15 less than the 
(early type) binaries considered here. 

As mentioned in Section [3] the estimated M p hot following 
from a comparison with SSP models, or from an IMF extrapola- 
tion from the number of RSGs as is done i n resolved clusters, is 
also affected by binarity toavies et alj 2009) The fraction of stars 
that is removed from our sample due to this effect is (1 — <?)/-. cor- 
responding to 45%, giving rise to Mj™ /Mphot ~ 2 for the values 
of Table [2] There is no reason, however, to expect that this would 
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preferentially affect clusters with low ratios Mphot/r^s and it can 
thus not cause the downward trend seen in Fig. [3] 

The values of the binary parameters used in the study (Ta- 
ble O are only indirectly based on observations since we have to 
correct the period distribution found for massive main-sequence 
stars to account for the reduced RSG phase of stars in tight bi- 
naries (section [2~TV Our assumption can be verified once the bi- 
nary fraction fso and the associated period distribution among a 
statistically significant sample of resolved SGs has been deter- 
mined. This could be done spectroscopically using a long time 
base (~few 100 — 1000 d). The recently discovered RSG clusters 
towards the Gala ctic centre jFiger et al1l2006l ; iDavies et al. I l2007l ; 
IClark et alj[200gh provide an excellent opportunity to do this. All 
three have approximately the same age as the extra-galactic clusters 
used here and their masses are relatively low (few times 10 4 M©) 
and have radii of a few pc, which according to our model places 
them in the regime where binaries dominate the measured velocity 
dispersion. The ratio Mffi /Mp hot was determined for two of them 
and is ~ 2 IDavies et alj|2008l) . lower than the extra-galactic clus- 
ters with comparable A/phot /V e ff (Fig. O, but still consistent with 
the lower 2a line of our prediction. This result is very sensitive to 
low number statistics since the number of RSG in these clusters is 
~20, so for fso = 0.25 we expect only a handful of binaries. 

iRitchie et alj J2009h present a spectroscopic multi-epoch sur- 
vey of luminous evolved stars in Westerlund 1. This cluster is 
slightly younger than the clusters considered here, thus its super- 
giants population is formed by more massive stars. They find a bi- 
nary fraction in excess of 40% among the 20 most luminous su- 
pergiants. Interestingly, they also find radial velocity changes of 
~15 — 25kms~ 1 in cool hypergiants due to photospheric pulsa- 
tions. Macro turbulence dispersions of 5 — 10 km s~ are also f ound 
for luminosity class II and III giants bv lGrav&Tone3 ( ll986h and 
ICarnev et al.1 j2008h . This is an additional complication in dynam- 
ical mass determinations of young star clusters containing massive 
giants. 

Our results are an important ingredient in the discussion on 
the importance of the early mass independent disruption, or 'infant 
mortality', of star clusters. The high velocity dispersions found for 
the clusters discussed here have been put forward as empirical ev- 
idence t hat many young (< 30 M yr) clusters are quickly dissolv- 
ing (e.g. lGoodwin & Bastian 2006, M08). We have provided argu- 
ments that the alleged super- virial state can largely be explained by 
orbital motions of binary stars. 

Early dissolution due to gas expulsion can still exist, but it 
probably occurs on much shorter time-scales (<< 10 Myr) than 
generally assumed. This idea is supported by the fact that the clus- 
ters considered here have densities of ~10 3 Mq pc -3 , correspond- 
ing to an internal crossing times of the stars of roughly 0.5 Myr. 
So these clusters have evolved for at least 20 crossing times. The 
crossing time in the embedded phase is much shorter than the cross- 
ing time at 10 Myr due to th e nonzero star fo rmation efficiency and 
the consequent expansion l lBastian et alj|2008h . The gas expulsion 
models show that clusters need about 20 initial crossing times to 
find a new virial equilibrium, or completely dissolve into the field 



(e.g. 


Goodwinlll997l;lGever & Burkerj200ll:lBaumaardt & Krouoal 


|2007 


). So at 10 Myr the super- virial state is undetectable and the 



clusters discussed here are therefore survivors of the gas expulsion, 
or 'infant mortality', phase. 
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